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Abstract 

We show that the l/p^ power corrections to the ultraviolet asymptotic solu- 
tions are allowed as consistent solutions of the coupled Schwinger-Dyson equa- 
tion for the gluon and (Faddeev-Popov) ghost propagators in Yang-Mills theory. 
This result supports the existence of the vacuum condensate (j4^) with mass 
dimension 2, as recently suggested by the operator product expansion and lat- 
tice simulations. We compare the solution with the result of operator product 
expansion. 
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1 Introduction 



The vacuum condensate is known to play the distinguished role for characterizing the 
non-perturbative aspect of the quantum chromodynamics (QCD). In fact, the gauge 
invariant vacuum condensates such as gluon condensate {-^flJ) of mass dimension 4 
and chiral condensate {ipip) of mass dimension 3 are well-known examples studied 
extensively so far. 

Recently, several groups have focused their attention to a novel type 

of vacuum condensate. In a series of lattice simulations 0, the gluon propagator 
and the running coupling constant of Yang-Mills theory (QCD without quarks) in 
the Landau gauge has been computed at large momenta, and it was shown that these 
data are compatible with the existence of the rather large 0(l/p^) correction to the 
perturbative predictions. In an operator product expansion (OPE) approach 
this correction is related to the existence of condensate in the Landau gauge.Q 
It was also claimed that this condensate might be related to instantons. The 
numerical investigations reached the sizable value for the vacuum condensate 
{s^^) ^ (L4GeV)2. 

In the Maximal Abelian (MA) gauge [0], on the other hand, other types of vacuum 
condensate of mass dimension 2, i.e., off-diagonal gluon condensate ^^^^ and off- 
diagonal ghost-antighost condensate C°'C°', have been proposed and investigated in 
relation to quark confinement These condensates can provide the non-zero mass 
for the off-diagonal gluons and ghosts, giving a natural explanation for the Abelian 
dominance confirmed by numerical simulations and quark confinement [§]. 

Although the composite operator ^/^^ is not a gauge invariant quantity, it was 
pointed out that its minimum along the gauge orbit can have a definite meaning [|l| . 
Moreover, it was shown 0| that is regarded as a Landau gauge limit of the more 
general composite operator of mass dimension 2 which is invariant under the Becchi- 
Rouet-Stora-Tyutin (BRST) transformation in the most general manifestly Lorentz 
covariant gauge. That is to say, the vacuum condensate of mass dimension 2 can be 
allowed to exist as a BRST invariant quantity. 

The main purpose of the present paper is to give another evidence of the exis- 
tence of 0(1/ p^) correction to the perturbative expressions which are expected to 
hold in the large momenta. In the present paper, we study the existence of such 
a power correction in the ultraviolet (UV) asymptotic region within the framework 
of the truncated Schwinger- Dyson (SD) equation for the gluon and ghost propaga- 
tors in Yang-Mills theory. In the previous paper [^], we have examined whether the 
coupled SD equation allows logarithmic corrections to the infrared (IR) asymptotic 
solutions with power momentum dependence predicted by Gribov and Zwanziger |jlO| 



and confirmed based on various methods; SD equation ||Tl],|12|,|13||lJ,[lg,[Ig], stochastic 



equation |T^, lattice simulations [T^. Such logarithmic corrections are proved to be 
absent in IR asymptotic solution to the truncated SD equation 0, in sharp contrast 
to the UV asymptotic solution. For the relationship between color confinement and 
the solution of the SD equation, see e.g., H]. 



^ In the previous work Q , the OPE was performed among others in the most general manifestly 
Lorentz covariant gauge which reduces to the Landau gauge in a special limit A — *■ in the notation 
of [0. Therefore, the OPE result for the gluon propagator in the Landau gauge Q was confirmed 
by Pi , while the OPE of the Faddeev-Popov ghost propagator was calculated anew in p] . 
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In order to show that both the gluon and ghost propagators can have the 0{l/p^) 
power correction to the UV asymptotic solution with logarithmic momentum depen- 
dence, we adopt the Ansatz for the gluon form factor F{p^) and ghost form factor 
GijP') defined by the gluon propagator in the Landau gauge and the ghost propagator 
with unbroken color symmetry in the SU{Nc) Yang-Mills theory, 

Dt.'ip) := 5^''PlAp)np')/p\ G^^ip) := (1) 
with A, B = 1, ■ ■ ■ , N^ — l and the transverse projection operator, P^^yip) ■= S^^ — ^^, 

N N' 

F{z) =Az^ Cn-z"" + E e"^^/^2^^4^\ (2a) 

n=0 i=\ 
N N' 

G{z) =Bz^ ^nz-"" + e-^"'"'zH^o\ (2b) 

n=0 1=1 

where0 cq = 1 = d^. Here we have introduced a variable ||19|,||, 



p2 

^ := u; In h z^^ (3) 

a 

where u; is a real number and a is the renormalization group invariant momentum 
scale defined by o := yU^exp[— 2 /^^ j^] for the renormalization scale /i using the /? 

function: P(g) := /^^^t^- The UV asymptotic limit — oo corresponds to 2; — > cx). 

2 Power corrections to asymptotic solutions of the 
coupled SD equation 




Figure 1: Diagrammatic representation of the truncated Schwinger-Dyson equations 
for gluon and ghost propagators. 



In this paper we deal with the renormalized version of the SD equation. 

^We can consider more general Ansatz which is equivalent to the sum of the logarith- 
mic and the power forms with possible logarithmic corrections: F{z) — Az'' X]!^=o '^"'^ " + 
'Y^=\ e~^^/'^z^* S^=o where 7^ := 7 + 7« and Kn^ := Aa'n Cn- The first correction factor 
in the above Ansatz corresponds to the higher-order logarithmic corrections in perturbation theory, 
while the second factor to the power corrections. 
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0) We introduce the renormalization constants for the gluon field Z3,_the ghost 
field Z3, the triple gluon vertex Z\ and the gluon-ghost-antighost vertex Z\ to write 
down the SD equation for the renormalized gluon and ghost form factors. 

First, we restrict our consideration to a version of the truncated SD equation 
subject to the following truncation, procedure and an approximation in order to 
write down the solvable SD equation. See Fig. |l] for the diagrammatic representation 
of the SD equation treated in the present paper. 

1) The full vertex function is replaced by the bare one. 

2) Contributions from the two-loop diagrams are neglected in the gluon equation, 
while all diagrams are included in the ghost equation. 

3) The gluon equation is contracted with the Brown-Pennington projection oper- 
ator RfjLvip) '■= ^pLv — 4:^^^, to remove the quadratic UV divergence coming from 



the tadpole term. 

4) The y-max approximation:^ F{{p — g)^) = F(max{p^, g^}), G{{p — qY) = 
(^(maxjp^, g^}), is adopted to avoid angular integration of the angle 6 where the 
angle 6 comes from the inner product p ■ q = ^/p^y/q^ cos 9 between the external 
momenta p and the internal (loop) momenta q which is to be integrated out. 

By taking the approximations 1), 2) and 3), the SD equation includes only two 
types of the bare vertex for the triple gluon and the gluon-ghost-antighost interactions. 
The non-trivial contribution in the gluon equation comes from only two diagrams, i.e., 
gluon loop and ghost loop, apart from the trivial tree gluon propagator. In particular, 
the projection 3) removes the tadpole diagram. On the other hand, the ghost equation 
includes a mixed loop composed of gluon and ghost. By the approximation 4), the 
coupled SD equations are reduced to the one-dimensional integral equations with 
integration variable q^ G [0, A^] where A is the UV cutoff. Possible improvements of 
the above approximations are discussed in the next section. 

Thus, the ghost SD equation reads 



G-\x) =Z,- ^AZi 



^^""^ rdyyG{y)+ f''^ ^F{y)G{y) 

Jo Jx V 



x"^ Jo " " ' Jx y 



(4) 



while the gluon SD equation reads 

^(^) r ^ 2^/ N 3G(x) r , , 1 f^' dy 



3 / dyy'Giy) + ^ / dyyG{y) + ^ [ ^G'{y) 
x'^ Jo /x^ JO 2 Jx y 



7F(x) r , , 17F(x) r , , 9F(x) r , , 

^ ^ / dyy^Fiy) - — / dyyFiy) - — ^ / dyFiy) 

Jo IX^ Jo QX Jo 



'- I A/1/2 pr,,^ _ im^iz. / _ v-"^ / 

2x3 

Jx y 8 Jx y^ 



(5) 



■^This approximation was called Landau- Abrikosov-Khalatnikov (LAK) approximation or 
Higashijima-Miransky (HM) |2^] approximation. In QED, the angular integration in the coupled SD 
equation in the bare vertex approximation was first performed by Kondo, Mino and Nakatani [2^ 
and subsequently by Bloch and Pennington p3| in a numerical way. In the pure Yang-Mills theory 
(QCD without quarks), the exact angular integration in the bare vertex approximation was per- 
formed in an analytical way for the pure power solution as the IR asymptotic solution jT^ w ithout 
relying on this approximation. Recently, a more sophisticated procedure was invented ]15, 16 for 



the solution to be consistent with the renormalization effect of the triple gluon interaction vertex 
which is necessary to reproduce the correct coefficient of /? function in the UV limit. This issue will 
be discussed in the next section in more detail. 
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where we have introduced A : = 



167r2 ' 



X 



and y := 



The two of the four renormahzation constants, and Z^, can be ehminated by 
subtracting the equation at x = cr. The ghost equation reads 



3N, 



\Zi 



F{x) 



x^ 



dyyG{y) 



[X 



<y) + 



" dy 

X y 



F{y)G{y) 



(6) 



while the gluon equation reads 



F-\x) =F-\a) + ^XZ, 



— (X 



3 



7F(x 



G{x) 
x"^ 

''G\y) 



dyy^G{y) + 



3G{x) 
2x2 



dyyG{y) 



X 2y 



2x3 



dyy^F{y) 



17F(x) 



2x2 



+ T r - (X - a) - 7 

8 Jx y^ 



dy r,2 



X y 



dyyF{y) 
F'\y) 



9F(x) 



5X Jo 



dyF{y) 
(7) 



The coupled SD equation obtained in this way agrees with the version adopted by 
Atkinson and Bloch . It is known that this solution exhibits qualitatively the same 
IR behavior as the solutions of the truncated SD equations using the dressed vertex 
function improved so as to be consistent with the Slavnov- Taylor identity. Moreover, 
it has been shown that the leading logarithmic behavior is reproduced in the UV 
asymptotic solution so that it is consistent with the perturbation theory at one-loop 
order. In this paper, we consider the sub-leading contribution and power-corrections 
in the UV asymptotic solution. 

In the actual calculations, it is more convenient to rewrite the coupled SD equation 
in terms of the new variable z := uj\n- + z^j and C '■= culn - + z„. Then the ghost 
equation reads 



G-\z)-G-\z^] 



3iV, 



'-XZi 



(x 



(t; z 



dC 



ujoo UJ 

dC 



FiOGiC) 



and the gluon equation reads 

F-\z)-F-\z^) 

G{z)e-^'/'^ 



3 



r ^e^a^GiO + -G{z)e~'^/^ f 

J —LOCO UJ 2 J — 



c^oo UJ 



GiO 



— X 



o"; z 
7 



^G\C) 



UJ 



-F{z)e'^'/- r ^e^i/^F{C) - —F{z)e~^'/'' f 
2 J —uioo UJ 2 J - 

8 J—U)OD UJ 8 J Z UJ 



UJ 



FiC) 



-(x^a-z^ z^)-7 I —F\C) 

UJ 



(9) 
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In order to find the solution, we substitute the Ansatz into the above equations 
(§) and (^), then perform exphcitly the integration over (, and finally comparing both 
sides of the resulting equations to match the coefficients and the powers on both sides. 

In the process of integration, we use the integration formula [E^ 



a 




<i+b)Y[l + b,-az] (a 7^0) 



"(^"('-i'-'/.Z,,... (10) 

where r[c, x] is the incomplete gamma function with the asymptotic expansion for 
large r[c, x] = x^~^e~~^ 1 + J2'^=i . In particular, for large z and 

a 7^ 0, we can use the formula. 



i+EM"'n(i-^) 



=1 



In order to use the UV or IR asymptotic solutions in the integrand, we perform the 
decomposition: /(f dyf{y) = J^2 dyf{y) + dyf{y) for the UV limit x ^ A^{-^ oo), 
or dyfiy) = J^dyf{y) + dyf{y) for the IR limit x —* e {e is introduced to 
avoid IR divergence). The asymptotic solution can not determine the constant part, 
/(f dyfiy), but it is a p-independent constant and does not affect the exponent. 

The first observation is that the Ansatz for asymptotic UV solution without power 
corrections (a^^ = and h^^^ = 0^ can not satisfy the coupled SD equation, even if 
we include the logarithmic corrections Cn,dn 7^ 0. This fact is made sure by explicit 
calculations. This is because the l/(p^)^ [i = 1, 2, 3) power corrections for x := p^ 

^ 1 are generated from the term with whole range of integration [0, A^] which 
appears after separating the integral for large x A^(— 00) as 



x'^ 



. r dyy'-'F{y) = ^ f dyy'-'F{y) + ^ / dyy'~'F{y), (12) 

JO X*^ JA2 X*^ Jo 



in the RHS of the gluon equation F~^{x). Due to the term — dyF{y), a 
solution without a leading power correction 0(1 /p"^) is not allowed in the truncated 
SD equation in question.| Thus, the gluon equation involves at least l/p^ ,1/p'^ ,1/p^ 
power corrections, while the ghost equation involves at least 1/p^ ones. The leading 
power correction originates from the gluon loop in the gluon equation. However, 
the leading power correction is also generated by other terms after integration, once 
the Ansatz contains the power correction part. Therefore, the actual procedure of 
matching the relevant terms is more complicated.^ Keeping these observations in 

^It is pedagogic to see whether it may happen that the integral f^ dyy^~^F{y) cancel the power 
dependence x^, i.e., dyy^^^F{y) = x^f{z) where /(z) is a polynomial of z. For this equation to 
be satisfied, F{x) must be of the form: F{x) = if{z) + Luf'{z), as obtained by differentiating both 
sides with respect to x. Therefore, if the solution F[x) is written as a polynomial of z in the whole 
range [0,a;] {x — > 00), the power dependence happens to cancel. This is impossible, since we know 
that the IR asymptotic solution does not have such a form, see e.g. ||]. 

^ In the fermion SD equation of QED and QCD, the similar situation arises if one tries to find 
the asymptotic solution, although the fermion SD equation can be linearized and be solved in the 
whole momentum region. However, we have checked that the asymptotic solution agrees with the 
correct solution T,{x) = Ce~^ z'^^^-^^^q RnZ~" with z := ln(x/AQ(^^), obtained by converting 
the fermion SD equation to the linear ordinary differential equation with appropriate boundary 
conditions |19|. 
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mind, we proceed the calculations. 

In the ghost equation, for the adopted Ansatz to be a solution, the following 
relations must be satisfied (see Appendix). 



-5 = 1+7 + 5, 



3iV, 



cu(l + 7 + 5)^ 



(5i = 7i - 1, 



.(1) 



9Nr 



(13a) 
(13b) 
(13c) 
(13d) 



where ( [13a|) and (|13b| ) for the leading logarithm come from the matching between 
the LHS and the last term of RHS, i.e., f F(C)G(C), while O and (^3^ for the 



(1) /.W 





power correction are due to all the terms in the RHS.0 (|13d| ) gives a ratio % / 

Similarly, the gluon equation implies that the parameters must satisfy the follow- 
ing relations. 



-7 



1 + 25 = 1 + 27, 



XZiAB' 



:_L_ + ^,Z,A'_1_, 

2cu(l + 25) 3 u;(l + 27)- 



-7 + 7i = 1 + 27 + 71, 



9/^ 



1 = !^XZ,A'- -. 

3 to'(l + 7 + 7i) 



(14a) 
(14b) 
(14c) 
(14d) 



Here (|14a|) and ( |14b|) for the leading logarithm come from the matching between 
the LHS and the last terms in RHS,[| i.e., ^.'^ fG\C) and ^F\C). The 

leading UV logarithmic behavior comes from the equal contributions of the ghost loop 
and the gluon loop. This is sharp contrast wit the IR power behavior whose dom- 
inant contribution comes from the ghost loop (ghost dominance) at least for the 
present truncation. While ( |14c| ) and ( |14dD for the power correction are due to a term 
-|e-^/'^F(2) Jl^^ f e^/^F(C) in the RHS.0 Therefore, m the gluon equation, the UV 

^ The relations ( |l3c| ) and ( |l3d[ ) are obtained as follows. The matching condition between LHS 
and RHS implies —6 + 61 = max{j + S + 7 + S + ji }. T he former choice —6 + 61 = j + S + 5i does 



not hold, since 7 + 2(5 = —1. Hence, the latter choice (13c) must be satisfied. Once this identification 
is made, (13d) immediately follows. 

^The relation (14a) is obtained from —7 = max{l + 2(5 (ghost loop), I + 27 (gluon loop)}. For the 
choice —7 1 + 2(5, 5 > 7 (ghost dominance), (14b) must be replaced by 1 = ^XZiAB^ 2^(7+215) • 
For this relation to be consistent with (|l3^), we have -6(1 + 2^) = (4/3) (1 + -f + 5) = -(4/3)(5 and 
hence S — —9/16. This leads to 7 1/8, contracting with (5 > 7. For the choice —7 = 1 + 27, (5 < 7 
(gluon dominance), i.e., ^ < 7 = —1/3 and (14b) must be replaced by 1 = ^"^^^^^ ^{1+2^) ■ "^^^^ 
choice is excluded easily , sin ce —7 = 1 + 26 from (|l3a|) leads to (5 = 7, contradicting with (5 < 7. 



Once the identification (14a) is made, (14b) immediately follows. 

® The relation (14c) is obtained as follows. By comparing both sides of the gluon equation, we 
obtain —7 + 71 = max{25 + (5i (ghost loop), 1 + 27 + 71, 7 (gluon loop)}. The two conditions from the 
gluon loop stem from the separation of the range of integration, — |e~^/'^F(z) Jl^^ P{C)i i^^to 

[zA, z\ and [— cjoo, za]- The matching condition is rewritten as 71 — max{7+25+(5i, I + 37+71, 27} = 
max{(5i — 1,71,27} — max{7i — 2,71,27} = max{7i,27} where (13a), (14a) and ( |l3cD were used. 
Therefore, the exponent 71 is determined from the gluon loop. When 71 > 27 = —2/3 — —0.666 • • • , 
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power corrections originate from the gluon loop contribution. We anticipate this re- 
sult, since the dimension two gluon composite operator can be inserted into the gluon 
loop only, see e.g., the second reference in ||^. 

We look for the UV asymptotic solution with leading power correction to the 
logarithmic solution deduced by perturbation theory. First, the mixed loop (|13aD 



and the ghost loop ( p.4a|) yield the same relation. 



7 + 25 + 1 = 0. 
However, the gluon loop ( |14a| ) implies 7 = —1/3, leading to 

5 = -1/3 = 7. 



(15) 



(16) 



This result is consistent with the relationship (3^) following from the comparison 
between exponents of the leading power correction. From ( p.3bD , uj is obtained as 

(9/4)A/'cAZiA52. (17) 



In order to reproduce the asymptotic freedom, oj must be positive, > 0. 
substituting this into ( 14b ), the renormalization factor Zi is expressed as 



Substitution of these relations into ( |14d[ ) yields the exponent 

7^ = -7 - 1 + 11/168 = -101/168 = -0.60119, 



By 



(18) 



(19) 



-2/3 



-0.666' 



in the footnote. Hence, 



which satisfied the condition 71 > 27 
another exponent is 

5i = 71 - 1 = -1.60119. (20) 

Finally, the ratio of two coefficients jcL'^^ is obtained from (|13d| ). It should be 
remarked that the absolute values of the coefficients for the leading power corrections 
can not be determined from the UV asymptotic solution alone, just as the coefficient 
of the logarithmic correction in the IR region can not be determined from the IR 
asymptotic solution alone 0. In order to determine the coefficients, we need to 
connect the UV asymptotic solution to the IR asymptotic solution and vice versa. 

The running coupling constant for the gluon-ghost-antighost interaction in the 
UV region is given by 



gl{v) ■.=g'Z,F{p')G\p) 



g^Z,AB^z"'+^^ [1 + (ci + 2di)z-^ + (^'^^^^ + b^^^ z^')e-'/'' + O ( 



-2 -4 

Z ,p 



(9/4)iV p 
log- 



167r2 



1 + (ci + 2di)z-^ + (4'^^^i + fei'^/^e-^/" + O {z-\p-^)' 

(21) 



( |l4cj ) is obtained from the integral in the range [za, z] where 1 + 7 + 71 > 0. When 71 < 27, we faU 
into a contradiction 71 = 27 to be excluded. When 71 ~ 2j, two integrals seem to give the same 
order of contribution. However, this case I+7 + 71 — 1 + 87 = leads to a new type of momentum 
dependence, i.e., double logarithm Inz ~ Inlnp^, which is beyond the Ansatz we adopted. Therefore 
the case 71 = 27 is excluded in the present analysis and should be investigated elsewhere. Once the 



identification (14c) is made, we obtain a'j^^ = ^XZiA^ 



9/8 
^1^(1+7+71)' 



, which implies (14d) 
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Another running coupling constant for the triple gluon interaction is defined by 



glip) ■.=g'Z,F\p') 

=g^ZiA^z^"' [l + 3ciz~^ + Sa'i^z^'^e-'/'^ + O 



(63/ll)iV, p 
log^ 



1 -1 



167r2 



1 + 3ciz-^ + 34'^;z^ie-"/" + O (z-^^p-" 



(22) 



This solution implies the UV asymptotic freedom. However, for the coefficient 
of the P function, it leads to Po = |iVc or Pq = ff^c which are different from 
the perturbative result Pq = yiVc. This disagreement was resolved by taking more 
involved renormalization prescription in p!T| , pr5[| . We improve the above analysis taking 
into account this modification. 



3 An improvement 

Within the presented class of truncation schemes, it is impossible to satisfy both the 
correct one-loop scaling and the Slavnov- Taylor identity. 

In order to ensure the correct one-loop scaling of the running coupling, Fischer, 
Alkofer and Reinhardt [15| proposed a two-parameter Ansatz, i.e., the gluon vertex 



renormalization constant Zi is replaced by the momentum dependent substitution 

Z,^Z,ip,q,ip-q)):= -^^^ F({p - q^^^ " (^3) 

This Ansatz is suggested using the scaling of the dressing functions extracted from 
the renormalization group equation |Tl] . 



In the UV limit, this Ansatz yields the correct one- loop behavior of the gluon 
loop, for arbitrary a, 6, as will be demonstrated below. The two parameters a and h 
must be chosen so that the momentum dependence of Zi is as weak as possible, since 
Zi would be independent of momenta in a full treatment of the coupled gluon-ghost 
system. 

On the other hand, the IR limit strongly depends on the choice of two parameters 
a, h. For 6 = 0, three different situation arise: For a < 0, the ghost loop is dominant 
in the IR limit and the gluon loop is sub-leading. For a = 0, the gluon loop reproduce 
the same power as the ghost loop in the IR limit. For a > 0, the gluon loop becomes 
the leading term in the IR limit. In this solution for the coupled gluon- 

ghost system has not been found. Explicit choices taken so far are as follows. The 
choice a = 6 = corresponds to the truncation scheme of |]T3. Another choice 



a = 36{< 0),b = is adopted in [^. For 6 7^ 0, a choice a = b = 36{< 0) is 



studied in the original work ||T5|. This choice is shown to minimize the momentum 
dependence of Zi, see Appendix A of [0. 

Now we adopt the Ansatz (^) to study the power correction in the UV region. 
Repeating the similar calculations as those in the previous case, we obtain another 
set of relations for the power and the coefficients. Note that the first set of relations 
from the ghost equation (|T^) is not affected by the substitution (^). The second set 
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from the gluon equation (113) is replaced by 



-7 = 1 + 25 (ghost loop) = 1 + 26 (gluon loop), 



1 = ^XZ^AB^ 



]_ I . . l-g-fe Ty2-a/S-b/5-2a-2b 

2o^(l + 25) "^3 coil + 26)'' 



-7 + 71 = 1 + 27 + 71 + 7(-2 -a-b) + 6{2- a/6 - h/6 -2a- 26), 



Nr. 



A 



9/8(-l -a- 5)yll-«-b52-a/5-6/5-2a-26 



3 cj[l + 7 + 7i + 7(-2 -a-h) + 6{2- a/6 - h/6 -2a- 26)] ' 



(24a) 
(24b) 
(24c) 
(24d) 



The relation ( P4a|) is the same as ( [13a| ). By substituting uo obtained from ( |13b|) 
into (|24bD , we obtain the ratio: 7/5 = -| [l - 14Zf^(A52+i/5) , By using the 
perturbative normalization A = B = 1, the ratio reduces to 7/5 = 26/9 which implies 

-9/44 = -0.204545. (25) 



7 = -13/22 = -0.590909, 6 
Returning to ( |13b| ) yields 

[ll/3)NrXZiAB'^ = {ll/3)NrX. 



(26) 



The relation (|24c|) is trivially satisfied under the relation ( P4a|) and does not give any 
restriction. For Oq' 7^ 0, the final relation (p4d|) with the substitution ( ^6|) gives the 
exponent of power correction: 

7i = -7 - 1 - 7(-2 - a - 6) - 5(2 - a/6 - h/6 -2a- 26) - 9/88(1 + a + 6), (27) 

whose value depends on the choice of a, 6: 

71 = -0.74845(a = 35 = 6), -0.323089(a = 35, 6 = 0), +0.102273(a = 0, 6 = 0). 

(28) 

^^Iq^" of the j3 function is correctly 



The result (|26| ) shows that the coefficient (3q 
reproduced by the identification /3ofi'^ = oj. 



4 Comparison with OPE result 

The UV asymptotic solution obtained above for the gluon and ghost propagators with 
power corrections should be compared with those obtained by the OPE. To one-loop 
order, the renormalized gluon propagator has the OPE [@,|25| 



Mp) + 



p2 A{N^ - 1) 



(29) 



where the Wilson coefficients are given by 



Aoip)=(9\p)/9'if^) 



In-^/ln. ^ 



13 

^ 22 



A, 



QCD 



A 



QCD, 



A2{p) =Nrg\p) (g\p)/g\fi) 



70 
"/30 



In -r-^ — / In ■ ^ 



35 
44 



A 



QCD 



A, 



QCD, 



(30) 
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with the coefficients of the (3 function and anomalous dimension of the gluon field 
and the composite gluon field: Po = ^^c, 7o = ^^c, = ff^c, 7o = lA^ — 
7o = jNc- As already mentioned above, the exponent of logarithmic dependence 
completely agrees with the perturbative result: 7 = 7o//3o- We find the exponent of 
logarithmic correction —(1 — ^) = —35/44 = —0.795455 from OPE. The best fit is 
obtained 7 + 71 = —0.913998 for the choice a = 35, 6 = among the cases treated 
in the previous section (The naive treatment yields the value —0.93452.) Thus, the 
logarithmic momentum dependence of the power correction obtained from the UV 
asymptotic solution of the SD equation can reproduce the prediction based on the 
Wilson coefficient A2 of OPE to one-loop order. 

On the other hand, the renormalized ghost propagator has the OPE ||^ 



Boip) + 



B2{p) {^Rm 

p2 4(iV2-l) 



(31) 



where 



Bo{p)={9\v)/9\^^ 



70 
/3o 



In-^/ln. ^ 



_9_ 

" 44 



A, 



QCD 



A 



QCD, 



2a 



In 



p 



-/In- 



_9_ 

22 



(32) 



7a2 



Aqcd ^qcd J 

— 70 = yAi'c. Hence, the solution 
-1.913998 from the OPE resuh 



with /?o = yAT,, 70 = f AT,, 7^2 = ffAT,, ^0 
of SD equation exhibits large deviation 6 + 61 
— (1 — = = —0.409091 for the ghost propagator. This is the issue to be 
resolved in the near future. 

The running coupling constant for the gluon-ghost-antighost interaction is given 

by 



glip) =g'Z,F{p')G'{p) 



g^ZiAB^z'+^^ [1 + (ci + 2di)z-^ + {a^^^z^^ + b^^ ^ z^')e-'/'' + O {z-^p-") 
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tip) [1 + (ci + 2di)z-^ + (4'^z^i + 6f,'^z^0e~"/" + O , (33) 

while the running coupling constant for the triple gluon interaction is modified as 

glip) =g^ZrF\p') = ^2^(p2)l-a-5^(^2)2-a/5-6/5-2a-2b 



_g2 j^l-a-b Q2-a/5-b/S-2a-2b ^■~l{l-a-b)+5(2-a/5-b/5-2a~2b) 

X "1 + ?,ciz~^ + 34^^;z^ie-"/^ + O (z~^,p 
--glerM [1 + Sci^-^ + ?>a^^h^'e-'l^ + O ( 



2 -4 

z ,p 



(34) 



where g^^rtip) ■= I 



(ll/3)Afc 



log- 



^QCD 



Two running coupling constants do not agree 



to each other beyond the one-loop expression, since the Slavnov- Taylor identity is 
sacrificed by the bare vertex approximation. 

The OPE for the coupling constant of triple gluon interaction shows 

g'ip) = glertip) {1 + ^ (Hp/^qcd))'^'''^''^^''-' {KK)/p'} , (35) 
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where (7o+7o)//5o-l = -9/44 = -0.204545 and i? := ^^Jf^^ Kqcd))"^'''^^'^'^^ 
For the running couphng, therefore, the agreement for the exponent 71, 5i between 
the SD equation with the OPE is not so good. This is because in our truncation 
of the SD equation the bare vertex approximation is adopted, while in the OPE the 
vertex correction was taken into account to the leading logarithm. 



5 Conclusion and discussion 

We have proposed a new Ansatz for obtaining the asymptotic solution of the cou- 
pled SD equation for the gluon and ghost form factors (or propagators) in Yang- 
Mills theory. We have found a set of consistent UV asymptotic solutions with the 
leading power corrections 0(l/p^) for the gluon and ghost form factors: F{z) = 

Az^ En=o Cn^z"" + e'^/'^z^+^i ttfj^^ G{z) = Bz^ En=o ^^n^"'' + e-^'"^ z^+H^q\ within the 
framework of the truncated SD equation of Yang-Mills theory in the Landau gauge. 

Within the present truncation of the SD equation, we have calculated the ex- 
ponents of logarithmic dependence in the power correction term 0(l/p^). The mo- 
mentum dependence of the power correction part in the gluon propagator agrees well 
with the OPE result by taking into account a suitable modification for the triple gluon 
vertex, while the ghost propagator does not show good agreement. In this paper we 
have adopted the bare vertex function and the simple modification for truncating the 
SD equation. Therefore, it will be important to see whether our qualitative result 
remains true even after the vertex function is improved so as to satisfy the Slavnov- 
Taylor identity. This might give a clue to resolve the discrepancy between the OPE 
and the SD equation especially for the ghost propagator. 

Thus, the SD equation suggests the existence of the leading power correction to the 
UV asymptotic solution. The existence of the leading power correction is consistent 
with the existence of vacuum condensate with mass dimension 2, i.e., (=2^^) 7^ 0, as 
argued recently by several groups It turns out that the existence of power 

corrections with additional logarithmic dependence of momenta leads to the existence 
of sub-leading terms with logarithmic dependence. 

A way to determine the numerical value of non-vanishing coefficient Oq"* or 69^^ and 
thereby the numerical value of the condensation (^^) is to perform the numerical 
calculation which enables us to connect the UV solution into the IR one. It will 
clarify how the UV solution with the vacuum condensate (^^) is connected to the IR 
solution in which the new mass scale Mqi e.g., the Gribov mass (scale) is expected 
to appear. Another way is to require the (Borel) summability of the infinite series 
N = 00 where the consistency of the theory in presence of IR renormalon can be 
used to estimate the numerical value of the vacuum condensate (^^). It is important 
to notice that the mixing between the leading power correction and the logarithmic 
parts does not occur. These issues will be studied in a subsequent paper. 

It is in principle possible to extend this work to obtain the sub-leading power 
correction to the order 0(l/p^), corresponding to the vacuum condensate of mass di- 
mension 4, e.g., {^^^) 7^ 0. Finally, it will be interesting to study the UV asymptotic 
solution in other covariant gauges, e.g., the Maximal Abelian (MA) gauge [0 where 
a different type of vacuum condensate of mass dimension 2 is expected to occur, i.e., 
simultaneous Bose-Einstein condensation of off-diagonal gluon and ghost 
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A Integration of SD equation 

After performing the momentum integration, the ghost equation reads 



M M 



N=0 



N=0 



+ 7 + 5 — n — m) 



1+7+25— n—m 



+ e 



M 



'- €=0 i=0 

^ ^=0 ^ ^ ^ i=0 

M £-1 

+ c„d„aW^^+^^+^-"-"^ E ^"'(-'^)' n + ^ - 7 - 5 - 71 + ^) 

€=0 i=0 

M £-1 -, . 

+ C„d^6W ^7+25+51 -n-m ^ ^-^(_^)^ + ^ _ ^ _ 5 _ + ,) l 

^=0 i=0 -I 

(36) 

where we have defined v := XZiAB^. The gluon equation reads 



M 



M 



Ar=o 



A?"=0 



n,m=0 



1/ E 



_ „l+7+2(5-n-m J j 



1 



1 



-n-m) 



_|_ g-z/a;^7+2(5+(5i -n-m 



^ ^=0 ^ ^ ^ i=0 
i X - I UJ 



M 



e-1 



+ dndn^bg^ - E n - - '^l + ^) 

^ ^=0 i=0 

^ e=o ^ ^ ^ i=o 



1 1 ^ 

(37a) 
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n,m=0 



37-n-m 



" £=0 ^ ^ i=0 ^ £=0 ^ ^ ^ i=0 

-1 



9 



M 



7 



M 



£=0 

e-1 



" £=0 1=0 ° e=o 



Z [—UJ 



yYl{m + n-2^ + i) 



i=Q 



+ Z 



1+37— n—m 



Ct7 



_|_ g-2;/i^^37+7i-n-m 



Oi)(l + 2^ — n — m) 



iE--^(l) ^(--7-7l + ^) 





17 
2 



Cr,,Crr,a:;' ~ E ^ ^o;*^ ]J (m - 7 - 7i + ' 







j=0 



17 

^ £=0 i=0 



,9 



M 



£-1 



- c„c^aW- E ^ n("^ - 7 + ' 

O £=0 i=0 



7 ( 

+ c„c^(aW + aW) ^ 2;-M -- + m - 27 - 71 + i) 



9 



_ _p-2/u)^l+37+7i -n-m 



£=0 

M £~\ 

CnCm{a^n^ + aW)7^ ^-^(-u;)^ + m - 27 - 71 + i) 
£=0 j=0 
1 1 



(1). 



AT, 



M 



n=0 



(37b) 



^^z/. The above integrations were 
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where we have defined v' := \Z\J^ = ^^i^ — 28 
performed under the assumption that 1 + 7 + 5, 1 + 27, 1 + 25, 1 + 7 + 71 are not 
integers. In the LHS, the coefficients are obtained as 

Go=l, Gi = -di, G2 = -d2 + dl, G3 ^ -ds + 2did2 - dl, 
G4 = — d4 + 2did3 — 3dld2 + di,- ■ ■ , 

Fo =1, Fi = -ci, F2 = -C2 + cl F3 = -C3 + 2ciC2 - cl 



C4 -|- 2C]^C3 — 3c^C2 ~h 



Q^o^^-b'o', Qh"^d,{2b\,'>-b[''), Pr^-a'o', Pr^c^{2a\;>-a['>), 



(1) 



,(1) p(i) 



,(1) 



(38) 
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